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Abstract 

A new type of self-similar hot viscous radiative accretion flow onto 
a rapidly spinning neutron star has recently been discovered. This 
"hot brake" flow forms in the two-temperature zone (close to a central 
object), but at a sufficiently low accretion rate and a high spin it may 
extend in the radial direction beyond ~ 300 Schwarzchild radii into a 
one-temperature zone. When the spin of the star is small enough, the 
flow transforms smoothly to an advection-dominated accretion flow. 

The properties of the hot brake flow are rather exceptional and 
surprising. All gas parameters (density, angular velocity, temperature, 
luminosity, angular momentum flux) except for the radial velocity are 
independent of the mass accretion rate; these quantities do depend 
sensitively on the spin of the neutron star. The gas angular momen- 
tum is transported outward under most conditions, hence the central 
star is nearly always spun-down. The luminosity of the hot brake flow 
arises from the rotational energy that is released as the star is braked 
by viscosity. The contribution from gravity is small, therefore the ra- 
diative efficiency may be arbitrarily large as M — > 0. We demonstrate 
that the flow is also convectively stable and is unlikely to produce a 
strong outflow or wind. 
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The hot brake flow is cooling-dominated (via Bremsstrahlung) and, 
hence, might be thermally unstable. The analysis of thermal con- 
duction in the hot gas shows that thermal transport is collisionless 
(non-Spitzer) and occurs via free streaming of electrons along tangled 
magnetic field lines. We find that conduction is strong enough to make 
the flow thermally stable. 

The very fact that the density, temperature and angular velocity of 
the gas at any radius in the hot brake flow are completely independent 
of the outer and inner (except for the star spin) boundary conditions 
implies that the flow cannot be smoothly matched to a general exter- 
nal medium as well as to general conditions on the star surface. We 
demonstrate that there are two extra self-similar solutions: one bridges 
the gap between the original solution and the external medium, and 
another represents a boundary layer between the bulk of the flow and 
the star surface, in which the gas temperature rapidly drops while the 
density builds up. 

Finally, we briefly discuss that a hot brake flow may form around 
other rapidly spinning compact objects: white dwarfs and black holes. 



1 Introduction 



Accretion flows around compact objects frequently radiate significant levels 
of hard X-rays, indicating the presence of hot optically-thin gas in these 
systems. This has motivated the study of hot accretion flows around compact 

stars. 

Zeldovich fc Shakura (| 69h and lAlme fc Wilsonl f)l973h sidered spher- 



ically free-falling plasma impinging on the surface of a neutron star (NS). 
They calculated the penetration depth of the falling protons and made pre- 
liminary estimates of the radiat ed spectrum. Their i d eas were followed up b y 
a n umber of later aut hors, e.g., Turolla et al. ( 1994 b Zampieri et al. ( 199ot . 
and Zane et al. ( 19981 ). who carried out more detailed computations of spec- 
tra. 

A fluid approach to s pherical accretion onto a NS was pioneered by 
Shapiro fc Salpeter ( 197ot . who worked out the structure of the standing 
shock in a spherical flow and computed the two-temperature structure of the 
post-shock gas and the resulting spectr um. The equival e nt pro blem for an 
accret ing white dwarf was analyz ed by Kvlafis Sz Lambl ( 1982| ). In related 
work, IChakrabarti fc Sahul ()1997l ) described the hydrodynamics of spherical 
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accretion onto black holes and NSs, but without including radiation pro- 
cesses. 

All of the studies listed above involve inviscid spherical flows and include 
shocks of some kind, caused by spherically accreting matter crashing on the 
surface of the accreting star. However, if the accretion flow has angular mo- 
mentum and viscous transport, there are strong reasons why a shock is not 
expected. A shock introduces a causal discontinuity between the accreting 
star and the inflowi ng gas. While such a discontinuity is not a prob lem 
for a spherical flow, Pringld ( 1977 . see also Popham &: Naravan 1992) ar- 
gued that it leads to a ser ious physical inconsistency if the flow is rotating. 
Popham fc Naravan (1992) showed that a causal viscosi ty prescription leads 



to con sistent viscous accretion solutions without shocks. IPopham fc Sunvaev 
(2001J) have calculated detailed boundary layer solutions for accretion onto 
a neutron star, but these solutions are based on the cool thin accretion disk 
model and are not relevant for the hot flows we discuss here. 

One situation in which a shock is possible with a rotating flow is if the 
surface of the accreting star lies inside the marginally s table orbit of a thin 
disk. One then expects "gap accretion" with a shock ((Kluzniak &: Wilson 



19911 1 . and there is no inconsistency in having a shock. The marginally stable 



orbit appears not to play an imp ort ant role in hot quasi-spherical flows (see 
Naravan. Kato. fc Honma 1997 and Chen. Abramowicz fc Lasota 1997^ 1 . It 



is, therefore, not clear that one would necessarily have a shock with a hot 
flo w even if the accreting star were ve ry compact. 

iDeufel. Dullemond fc Spruitl (|200lh modified the model of lZeldovich Sz Shakural 
(1969) by considering a rotating advection-dominated accretion flow (see be- 
low) around a NS. Their model represents an improvement on the earlier 
work since it includes angular momentum and viscosity, but it still invokes 
a shock of some kind, since there is a discontinuity at the radius where the 
hot ADAF meets the surface of the NS. 



Fo llowing the seminal work on two-temperature accretion flows bv lShapiro. Lightman. fc Eardlc 
( 1976L the SLE solution), other hot solutions were discovered to describe ac- 
cretion onto black holes: the advection-dominated accreti o n flow (ADAF) 
(llchimanJll977tlR.ees et alll982tlNa"ravan fc Yill99lll995allbl:lAbramowicz et al. 
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, , . the advection-dominated in flow-outflow solution (ADIOS) ([Blandford &: Begelman 

,19991 : iNaravan fc Yilll994 Il995ah and the convection-dominated accretion 
flow fCDAF) (|Naravan et alitor iQuataert fc Qruzinovll200nh . All of these 
solutions describe rotating flows with viscosity and angular momentum trans- 
fer. The relevance of the solutions for accretion onto a NS is, however, un- 
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clear. 

Recently, iMedvedev fc Naravanl (|2001al ) discovered a rotating shock-free 
solution of the viscous fluid equations that corresponds to hot quasi-spherical 
accretion onto a rapidly spinning NS. We refer to this solution as a "hot 
brake flow" since the gas is hot (the temperature is nearly virial) and it 
spins the NS down (in previous works it was called the "hot settling flow"). 
The solution could equally well be described as a "hot atmosphere" since 
the solution is to first approximation static, and accretion represents only 
a small perturbation on the static solution (as is probably true for any set- 
tling flow). To our knowledge, the hot brake flow is the only solution for 
accretion onto a NS presently available that does not involve a discontinu- 
ity near the surface of the star. The hot brake flow should not be confused 
with a boundary layer which for ms in the very vicinity of the stellar surface 
(e.g., Naravan Sz Popham 19931 ) and which is characterized by a high gas 
density and steep spatial gradients of physical parameters. The hot brake 
flow forms above the boundary layer and extends radially to a large dis- 
tan ce, typically thousands of stellar radii or more . Following up the work 
by (jDavies fc Pringlelll98lh . Ilkhsanovl (|200lL l200ah has recently presented a 
subsonic hot accretion flow around a magnetized neutron star in the propeller 
state. The main difference between the two cases is that in the hot brake 
flow, heating and cooling balance each other, whereas in the subsonic pro- 
peller the heating rate of the accreting gas due to viscous dissipation is much 
larger than the radiative cooling rate. Hence, the latter solution shows some 
resemblance to the advection-dominated (or convection-dominated) flow. 

The hot brake flow exists at rather low accretion rates, smaller than a 
few percent of Eddington. The flow is subsonic everywhere (which is why it 
does not form a shock near the NS surface). Because the accreting gas has a 
low density and high temperature, the particle mean free path is larger than 
the local radius of the flow and the gas is essentially collisionless. Viscosity 
plays a very important role; indeed, the flow is powered by the rotational 
energy of the central accretor which is braked by viscous torques. A very 
interesting property of the flow is that, except for the inflow velocity, all 
gas properties, such as density, temperature, angular velocity, luminosity, 
and angular momentum flux, are independent of the mass accretion rate, as 
might be expected from the earlier comment that accretion behaves like a 
minor perturbation on an i ntrinsically static solution; the flow properties do 



depend on the NS spin fsee IMedvedev &: Naravanll2001al for more details). 



In this paper we present the analytical self-similar solution describing 
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this hot brake flow, study its stability and observable properties. We also 
present the solution that matches the hot break flow to an arbitrary external 
medium and the boundary layer solution that matches the flow to the star 
surface. All our theoretical results are confirmed with numerical solutions of 
the appropriate set of hydrodynamic equations. 



2 Basic considerations 
2.1 The mathematical model 

We consider gas accreting viscously onto a compact spinning object with a 
surface. The central object has a radius i?*, a mass M* = mMg un , and an an- 
gular velocity £7* = s£l K (R*), where ^Ir(R) = (GM^/R 3 ) 1 / 2 is the Keplerian 
angular velocity at radius R. We measure the accretion rate in Eddington 
units, rh = M / M^di and the radius in Schwarzchild units, r = R/R g , where 
^Edd = 1-39 x 10 18 m g s _1 (corresponding to a radiative efficiency of 10%) 
and R g = 2GM*/c 2 . We use the height-integrated form of the viscous hy- 
drodynamic equations ( Ichiman] 119771 lAbramowicz. et aL 19881 Paczvriskil 



199lUNaravan Yilh 994ft : 



M = 4nR 2 pv, (1) 

j-MVLR 2 , (3) 

P vT v~TB = — 7T To - vc 6 pJ^ = 9 - <?Coui, (4) 

P* vT *TB = — -~TB ~ vc Qe-jft = <?Coui - q , (5) 

where p is the mass density of the accreting gas, v is the radial infall ve- 
locity, Q is the angular velocity, c 2 = c 2 {6 p + 6 e m e /m p ) is the square of the 
isothermal sound speed, T p e are the temperatures of protons and electrons, 
6 P ,e = kBTp^/m Pfi c 2 are the corresponding dimensionless temperatures, v is 
the gas viscosity, J is the rate of accretion of angular momentum, s p and s e 
are the specific entropies of the proton and electron fluids, p e ~ {m e /m p )p is 
the mass density of the electron fluid, 7 P and 7 e are the adiabatic indexes of 









4 ™" B dR = 
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protons and electrons (which, in general, may be functions of T p and T e ) , and 
q + , q~, and gcoui are the viscous heating rate, radiative cooling rate, and en- 
ergy transfer rate from protons to electrons via Coulomb collisions, per unit 
mass. We have assumed that all viscous heat goes into the proton compo- 
nent. Equations (HJ-© describe the conservation of mass, radial momentum, 
angular momentum, proton energy and electron energy, respectively. 

We employ the usual a prescription for the kinematic coefficient of vis- 
cosity, which we write as 

v = ac s H f& ac s R. (6) 

Often, in accretion problems, one makes use of the relation H = c s /Qk and 
writes v = ac 2 /flK - This prescription is equival ent to u = ac s H rs ac s R 
in the regime of the iMedvedev fc Naravan ( 2001a ) hot brake flow. However, 



in the outer regions of the flow, where the other solutions described in the 
following sections appear, H is much less than c s /Qk, and v = oic 2 /VLk is 
not a good approximation. Equation (JHJ) is a superior prescription and is 
physically better motivated over a wide range of conditions (so long as the 
flow is quasi-spherical). 

We assumed that the flow is hot and quasi-spherical, which generally 
requires a low mass accretion rate (see Naravan. et al.lfl997h . The accreting 



gas has nearly the virial temperature, i.e., c 2 ~ GM^/R ~ (Q K R) 2 , and 
the local vertical scale height H = c s /Qk is comparable to the local radius 
R. We may then use the height integrated hydrodynamic equations for a 
steady, rotating, axisymmetri c flow, and for simplicity we may set H = R 



(jMedvedev fc Naravanll2001afl . 

In the case of accretion onto a NS we expect the flow to slow down as it 
settles on the stellar surface, and we expect the density in this settling zone 
to be significantly higher than for a black hole (BH) accretion. The increased 
density would cause more efficient transfer of energy from protons to electrons 
via Coulomb collisions and more efficient radiation from the electrons. As 
we show below, this leads to a flow in which q + , q~ and gcoui are all of the 
same order, which is very different from the case of a BH ADAF, where 
q + ^> q~ , gcoui- Another feature of the settling zone, again the result of the 
large density, is that optically thin bremsstrahlung cooling (which is sensitive 
to p) dominates over self-absorbed synchrotron cooling. We therefore neglect 
synchrotron emission in our analysis. For simplicity, we neglect also thermal 
conduction (we will include this effect in the discussion of thermal stability). 
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The temperature of the gas determines the efficiency of the Coulomb en- 
ergy transfer from the protons to the electrons and the rate of Bremsstrahlung 
cooling of the electrons. The balance between them defines whether the gas 
is in the two temperature regime (when the temperatures are high and the 
Coulomb collisions are very rare) or in the one-temperature regime (then 
temperatures are lower). 

The two-temperature regime occurs closer to the central object, where the 
virial temperature high enough and the electrons become relativistic. In the 
the inner region of the flow, < R < 100.R*. we expect a two-temperature 
plasma, with T p > T e , in which the electrons are relativistic and the protons 
are non- relativistic: 9 e >> 1, 9 P <C 1. The viscous heating rate of the gas, the 
energy transfer rate from the protons to the electrons via Coulomb collisions, 
and the cooling rate of the electrons via bremsstrahlung emission are given 
by 

g+ = ^W^V, (7) 



d R , 

TTt 

<?Coui = Qcouip 2 /, Qcoui = 4vrrglnA— e — , (8) 
Op mi 



v 

Q S ,RP%, Qff,R = 48«/r^> (9) 

m p 



where Oif is the fine structure constant, r e is the classical electron radius, 
In A ~ 20 is the Coulomb logarithm, c 2 s ~ c 2 9 p , and we have neglected 
logarithmic corrections to the relativistic free-free emissivity. The subscript 
"R" in Qff,R denotes relativistic Bremsstrahlung. 

For R > 100-R*, both protons and electrons are non-relativistic and have 
nearly the same temperature T p — T e T p , T e . The free-free cooling takes 
the form 

q- = Q S ,NR P 2 el /2 , Q ff , N R = 5V2n-^a f a T ^, (10) 

m A p 



where <jt is the Thompson cross-section, and the subscript NR stands for 
non-relativistic. Since the gas is effectively one-temperature, we have 

q + ~ q~ (11) 

because the Coulomb transfer rate, gcoui, is proportional to (T p — T e ) and can 
adjust itself to have the right magnitude with small adjustments of the two 
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temperatures. In the non-relativistic regime, 9 e < 1, 9 P ~ (m e /m p )9 e <C 1, 
and the Coulomb transfer rate is 

3 m e a T c 2 kT p -kT e 

fc - , = 7s^F lnAp ^v^ (12) 

From the condition gcoui — qT it follows that 

kT p - kT e ~ I^^-yTFv^c^^- 1 ^, (13) 

which is exponentially small for 9 e < 1. 

The set of equations (JU)-© must satisfy certain boundary conditions at 
the neutron star. First, as the flow approaches the surface of the star at 
R = i?*, the radial velocity must become very much smaller than the local 
free-fall velocity. Second, the angular velocity must approach the angular 
velocity of the star 

The radius of the star, and its spin, are the two principal boundary con- 
ditions applied at the inner edge of the accretion flow. We assume that the 
star is unmagnetized, so there are no magnetospheric effects to consider. Two 
outer boundary conditions, namely the temperature and angular velocity of 
the gas, are determined by the properties of the gas as it is introduced into 
the accretion flow on the outside (we will discuss this in subsequent sections). 
An additional important boundary condition is the mass accretion rate M, 
which is determined by external conditions and which we take to be constant. 

2.2 Numerical solution 

The system of equations given above with appropriate boundary conditions 
has been solved using the relaxation method on a highly non-uniform grid 
(in order to resolve the thin boundary layer, where the density rises by few 
orders of magnitude). 

We employ the gravitational potential of iPaczvhski fc Wiita f)l98nh to 



mimic the effect of strong gravity near the NS surface. In this potential the 
Keplerian angular velocity takes the form 



- GM * (14) 

Note that the analytical work presented in the following sections is based on 
a Newtonian potential. 
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We specify the boundary conditions as follows. We take the outer bound- 
ary of the flow to be at R out = 10 6 R g . At this radius we specify that the 
angular veloci t y is e qual to its value in the self-similar ADAF solution of 



Naravan fc Yil (J199J), and that the proton and electron temperatures are 
both equal to the self-similar ADAF temperature. We assume that the ac- 
creting star is a lMs un neutron star with a radius R* = 3R g = 8.85 km, 
unless stated otherwise. At R = R*, we specify the value of the NS spin 
parameter, s = Q*/Qk{R*), and we require the proton temperature of the 
flow to be T = few x 10 7 K -C T viria i. (We do not assume that the electron 
and proton temperatures are equal, but in fact they are equal.) We do not 
constrain the density of the gas in any way at either boundary. 

The numerical problem as posed here has a family of solutions charac- 
terized by three dimensionless parameters: the mass accretion rate rh (in 
Eddington units), the NS spin s (in units of the Keplerian angular velocity 
at the NS surface), and the viscosity parameter a. The angular momen- 
tum flux J, or equivalently the dimensionless flux j = J /MQ K (R*)Rl, is an 
eigenvalue of the problem. 

It is known that a hot flow with low rh around a black hole is an "advection- 
dominated" accretion flow (ADAF). It has a two-temperature structure for 
R < 300R g and is very hot (nearly virial) for all R. This solution has indeed 
been obtained numerically for rh < 0.01 and s < 0.01. An interesting fea- 
ture of the s = 0.01 ADAF-type solution is that it consists of two distinct 
segments. For large radii (beyond R ~ 20R g in Fig. the flow corresponds 
to the standard ADAF discussed in the literature, with the scalings 

p oc r~ 3/2 , c 2 s oc r" 1 , Q oc r~ 3/2 , v oc r' 1/2 . (15) 

However, at smaller radii, the numerical solution indicates the presence of a 
second advection-dominated zone, a "s ettling ADAF " whic h was first seen 



in numerical calculations described in iNaravan &: Yil (J1994). This settling 
ADAF is seen in Fig. as a zone that lies between the boundary layer region 
and the outer standard ADAF, with different slopes for p and v. The radial 
extent of the settling ADAF zone may be quite large and, in general, depends 
on the flow parameters and boundary conditions. 

A self-similar model of the settling ADAF may be readily obtained as 
follows. In an ADAF, energy is not radiated, therefore q~ = 0. Close to the 
star Q ~ constant, therefore q + = 0. Equations (jlj), (0) then simplify to 
the condition of entropy conservation, ds/dR = 0, which yields c 2 s oc p 7_1 . 
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Figure 1: Profiles of density p (g cm -3 ), proton temperature T p (°K), electron 
temperature T e (°K), angular velocity f2 (in units of the Keplerian angular 
velocity at the NS radius i?*), and radial velocity v (in units of c) for accretion 
flows with a = 0.1, m = 0.01, 7 = 4/3 and s = 0.3 (solid curve) and s = 0.01 
(dotted curve). The self-similar slopes for an ADAF flow and a settling 
flow are shown for comparison. The long-dashed curves represent the same 
solution as the solid curve, but with ten times higher temperature at R ou t- 



10 



As the material settles on the star, its radial velocity decreases and we have 
v -C Vff, Q <C Q K . Then, from equation (J2J), it follows that the temperature 
of the gas is nearly virial. Other quantities are determined straightforwardly, 
so that we have 

1 2-1 2 7 -3 

p oc n- 1 , c s oc r , il ~ const., t> oc r t- 1 . (16) 

The infall velocity decreases with radius if 7 < 1.5, and increases if 7 > 1.5. 
To highlight the difference between the standard ADAF and the settling 
ADAF, we have chosen 7 = 4/3 in the solutions shown in Fig. ^ 



3 Hot brake flow 
3.1 Numerical discovery 

By solving the equations numerically for different values of s, we have found 
that the with increasing s the ADAF solution continuously transforms into 
a solution of another type For relatively rapidly rotating stars, with s > 0.1, 
the new hot break solution (which was referred in the previous works as the 
settling solution) is already well-established. The transition is not sharp, 
so it is difficult to identify a specific transition point s = s t at which the 
transformation occurs. Numerical experiments indicate that the value of s t 
(however it is defined) is not very sensitive to R ou t, 7, and m and is, roughly, 
s t ~ 0.04 - 0.06. 

The change of the nature of the flow as s is varied is illustrated in Fig. 
[0 The solid and dotted curves correspond to two solutions with s = 0.3 and 
s = 0.01, respectively, with all other boundary conditions being the same. 
We see that the solutions are markedly different from each other. This is 
most clearly seen in the profiles of density, where the s = 0.3 model has a 
logarithmic slope of —2, as appropriate for the cooling-dominated settling 
solution described in this paper, and the s = 0. 01 model has a slope of —3/2 , 



as expected for a standard self-similar ADAF ([Naravan fc Yilll994l Il995al ). 
There is a similar difference also in the profiles of the radial velocity, where 
the two solutions have logarithmic slopes of —1/2 and 0, respectively. 

Figure El shows representative solutions for a = 0.1 and a range of values 
of fa and s. The solutions clearly have three radial zones. For R > 300R g , 
there is a one-temperature zone in which the gas properties vary roughly 
as power- laws of the radius. For R < 300i? 9 , there is a second power- law 
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zone with a two-temperature structure. Finally, close to the NS, the flow 
has a boundary layer region. In this final region, the gas experiences run- 
away cooling, the velocity falls precipitously, and the density increases very 
rapidly. 

The solutions are suspect also in the inner region of the two-temperature 
power-law zone and in the boundary layer, where Comptonization is likely 
to be important. Outside these regions, however, the numerical solution is 
expected to be accurate. 

We note here some qualitative results: (i) the boundary layer always 
forms near the star surface, (ii) the transition radius where the boundary 
layer meets the outer settling flow is usually at R tr ~ -R*, (iii) the value of R tr 
depends on the NS spin. These results are i n agreement with the studies b y 
Titarchuk. Lapidus Sz Muslimov ( 19981 ) and Titarchuk fc Osherovich (1999), 



and are important for the interpretation of kHz quasi-periodic oscillations. 

An interesting property of the new solution is its insensitivity to the value 
of accretion rate (provided it is low enough). In our numerical solutions, 
curves corresponding to a given value of s and different values of rh coincide 
with one other to very good accuracy. This is best seen in the profiles of p 
and Q. Changing s causes an up/down shift of the curves but does not affect 
the slopes of the curves. The temperature profiles are sensitive to the spin 
s, especially for large values of s. The radial velocity varies approximately 
as v oc rh and is roughly consistent with v oc r° for s > 0.1. 

3.2 Self-similar solution 

We now find the hot break solution analytically. The hot brake flow forms at 
low mass accretion rates. Therefore, for simplicity, we set m = and omit the 
continuity equation. Thus, the gas configuration corresponds to a radially 
static "atmosphere." The motivation for this approximation follows from 
the observation that the density p, temper a ture T and the angular velocity 



of the gas in the iMedvedev fc Naravanl (j2001af ) self-similar solution are 



completely independent of rh. Only the radial velocity v depends on rh, and 
it is given trivially by the spherical continuity equation 

M , . 

(17) 



AttR 2 P 

Because of this, we do not lose any generality by setting rh = v = in 
the analysis; we may always introduce a finite rh and finite v after the fact. 
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Log R Log R 

Figure 2: Profiles of density p (g cur 3 ), proton temperature T p (°K), electron 
temperature T e (°K), angular velocity Q (in units of the Keplerian angular 
velocity at the NS radius Rns), an d radial velocity v (in units of c) for 
accretion flows with a = 0.1 and (m, s) = (0.01,0.3) - solid line, (0.0001, 0.3) 
- short-dashed line, (0.01,0.1) - medium-dashed line, (0.01,0.7) - long-dashed 
line. 
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When d ^ 0, the energy equation has an extra term corresponding to the 
advection of energy. In advection-dominated accretion flows, for i nstance, 



this term dominates over the cooling term q~ ([Naravan. et al.lll9 97^. In the 
present case, however, we consider a situation in which the advection term 
is negligible (which corresponds to low fa). 

We assume that the flow is highly sub-Keplerian, Q(R) <C Q K (R), so that 
the centrifugal support is negligible compared to the pressure support. The 
radial momentum equation then takes the following simple form, 

GM* 1 d{pc 2 s ) 



R 2 p dR 



;is) 



where we have used the fact that v ~ an d written the pressure as y = 
pc 2 s where c s is the isothermal sound speed. iMedvedev fc Naravanl (J2001a|) 
present a more complete analysis in which they do not assume that the 
rotation is slow. They then obtain an extra factor of (1 — s 2 ) in their equation, 
which propagates through to all the results. Since we ignore the factor, our 
analysis corresponds to the case of a slowly-spinning star: s 2 C 1. This 
approximation is made only to simplify the analysis, and all the results may 
be generalized for arbitra ry s. 



From the analysis by IMedvedev fc Naravanl (|2001af ) , we know that the 



accreting gas in our problem acts as a brake on the central spinning star and 
transports angular momentum outward through the action of viscosity. We 
therefore write the angular momentum conservation equation for the gas as 
follows, 

J = invpR 4 —— = constant, (19) 
dR 

where J is the outward angular momentum flux, and v is the kinematic 
coefficient of viscosity. This equation is exactly valid in steady state if fa = 0. 
When fa is non-zero, there is an additional term, MQR 2 , due to the flux of 
angular momentum carried in by the accreting gas. The key feature of the 
hot brake solution is that the latter flux is negligible compared to the outward 
flux from the star. Equation (|19|) is, therefore, valid even when fa ^ 0, so 
long as fa is small enough for the term MQR 2 to be negligible. 

We first consider the one-temperature regime. The energy equations for 
the electrons and protons reduce in this case to 

q + ~ q~. (20) 
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The self-similar solution for the hot brake flow is then straightforwardly ob- 
tained. The gas parameters have the following radial dependences: 



p = Pir - 2 , T = T ir -\ Q = Q ir - 3/2 . (21) 

The subscript "1" in the coefficients is to indicate that this is the first (hot 
brake) solution, to distinguish it from other solutions described later. By 
substituting the above solution in equations (JTBJ), (fT^j) and (J2TJ|) . we see that 
it satisfies the basic conservation laws. We may also solve for the numerical 
constants: 



2 n / \ 1/2 

as 9 / m, N 



* = c3Qff ' NR ' (22) 

W = ^f, (23) 

We note that if rh ^ then the flow has a small constant radial velocity: 

v oc r°, (25) 

as follows from equation (fT7j) . 

The angular momentum flux in the solution is given by 

j = -a 2 s 3 R 2 a K^(-) V2 (26) 

By assumption, this flux is much greater than the angular momentum flux 
due to accretion, which sets an upper limit on the m ass accretion rate for 
the solution to be valid ( Medvedev fc Naravan 2001a| ). 
The pressure in the flow is given by 

P = pc 2 s = Picli r ~ 3 = Pi r ~ 3 , (27) 

where c 2 sl = 2kT\jm v . If the accretion flow is immersed into an interstellar 
medium with some external gas pressure p ext , then the above self-similar 
solution describes the flow at radii r (pi/Pcxt) 1 ^ 3 , where the pressure 

P > Pext- 



15 



This solution can be generalized to the two-temperature regime. In this 
case the electrons and proton energies are governed by two separate equa- 
tions, which under our assumptions reduce to 

q + = gcoui and g C oui = q~, (28) 

respectively. The solution given by equations (J2TJ) remains unchanged, with 
T being the temperature of the proton component, 

T p = T pl r~\ (29) 
T pl = 271, (30) 

and the temperature of the electrons is 

T e = T el r- 1 ' 2 , (31) 

kT el = (9^^kTX /2 . (32) 
V Qs,R Trip J 

The obtained hot brake solution has the remarkable property that all the 
quantities are uniquely determined by a single parameter s — the dimension- 
less spin of the central object — specified on the inner boundary. The fact 
that the solution does not depend on the outer boundary condition in any 
way means that there is no simple way to match it to the external medium. 
Clearly, there has to be a second solution to bridge the gap between this 
solution and the external medium. We derive the bridging solution now. 



4 Matching the hot brake flow with external 
medium 

We assume that the spinning star is immersed in a uniform external medium 
with a density p ext , temperature T ext and pressure p ext . We seek an accretion 
flow solution that extends from the spinning star on the inside to the external 
medium on the outside. As we have seen above, there is the self-similar hot 
brake solution extending from the boundary layer R ~ outward through 
a large distance, at least a few hundred R* or more. However, this self- 
similar solution has the surprising property that the density, temperature 
and angular velocity of the gas at any radius are completely independent of 
the outer boundary conditions. Hence, the solution cannot be matched to a 
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general external medium. We resolve this paradoxical situation by showing 
that there is a second self-similar solution whic h bridges the gap between th e 
original solution and the external medium (see Naravan Sz Medvedev 2003). 



This new solution has an extra degree of freedom which permits it to match 
general outer boundary conditions. 

We consider next the gas that lies just outside the region of validity of 
the first self-similar solution described above. In this zone, the pressure is 
expected to be approximately equal to the external pressure p ext ' 

pc 2 s = pext = constant. (33) 

This condition replaces the hydrostatic equilibrium equation (fTSJ). while equa- 
tions (fT^j) and (ffilj) continue to be valid. In this region, we find that there is 
a second self-similar solution of the form 



p = P2T 



- 7 '\ T = T 2 r 7 '\ n = n 2 r- g /\ (34) 



where the label "2" refers to the fact that this is our second solution. 

To match the second and first solutions, we require that the fluxes of 
angular momentum in the two solutions must be equal; this yields the con- 
straint (3/2)p 1 fi 1 T 1 1/2 = (9/4)p 2 ^ 2 T 2 1/2 . Making use of this and the other 
equations, we solve for the numerical coefficients in equation (jMj) : 

a 3/2 s 3 35/2 (m X3/4 ^ 



92 P^f2^Uj Ql%> (35) 

a 1 / 4 ^/ 2 , , fm \ 1/8 r 2 
^ = 4^- (37) 



l/4p5/4 \ m I l/4 

The pressure in this solution is constant and equal to the external pressure, 
p ext , and the angular momentum flux is also constant and is equal to J in 
equation (j2SJ). If the flow has a small but nonzero accretion rate, m ^ 0, 
then its radial velocity varies as [see eq. fTTjl ] 

v oc r 3/2 . (38) 

Whereas the original hot brake self-similar solution has a unique profile for 
a given choice of s, we see that the second solution derived here has an extra 
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degree of freedom, namely the external pressure p ext . This extra degree of 
freedom solves the problem discussed above. Thus, th e full solution consists 



of two zones: an inner zone described by the first ([Medvedev fc Naravan 



2001aJ) solution (|2Tj) and an outer zone described by the second solution 
The radius r matc h at which the two solutions match is obtained by 



n/ l/3„2/3 ql/3 f m \ 1/6 



equating the pressures: 



r mat ch - p l/3 R l/3 2 7/6 [ m J ^s,nr- ( 39 ) 

The second solution matches the external medium at the radius r ext at 
which its temperature matches that of the medium. This gives 



a 3/7 s 6/7 3 5/7 / m \ 3 / 14 c 8/7 



^ vTdkT^RT 25/14 UJ m p Q 3 s %' m 

If we wish we could also write this in terms of the external density by making 
the substitution fcT ext = m p p ext /2p ext . 



5 External medium solution 

For completeness, we presen t here the solution inside the external medium 



ror c ompletene ss, we presen t nere tne solution msiae tne external meaium 
()Naravan fc Medvedevl liffll . By assumption, the external medium has a 



uniform temperature and density, and a uniform rate of cooling. To main- 
tain equilibrium, there has to be some constant source of heat that exactly 
compensates for the cooling. We assume that such a source of heat exists 
(e.g., cosmic rays). The rotation Q is non-zero, but it decays rapidly outward. 
The small amount of rotation helps to transport the angular momentum flux 
from the star out into the external medium. Solving the angular momentum 
conservation law (3), we obtain the following solution 

p = p ext , T = T cxt , Q = Q 3 r~ 4 , (41) 

where 

o5/2 1/2 5 

^3 = a , 3 (^ ext )" 1/2 , (42) 

and p ext = 2kT ext p ext /m p . For m^O, the velocity scales as r~ 2 . 
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We confirm the existence and the structure of our matching solution (and 
the ex ternal medium) using the same numerical model ((Naravan fc Medvedev 
2003). In the calculations, the flow was taken to extend from an inner ra- 



dius -Ri n = 3 R g to i?out = 10 7 R g . The mass accretion rate was taken 
to be low, rh = 2 x 10 -5 , in order that the flow should correspond to the 
regime of the hot settling flow solution. We took the viscosity parameter 
to be a = 0.1 and set the spin of the star to be s = 0.3 (i.e., 30% of the 
Keplerian rotation at the stellar surface). We took the other inner boundary 
conditions to be the same as in MN01. At the outer boundary, we speci- 
fied the temperature and density of the external medium. Figure El shows 
four solutions. The external temperature is kept fixed at T(R ext ) = 10 8 K 
in all the solutions, but the external density varies by a decade and a half: 
p(R ext ) = 2.5 x 10 9 , 8.1 x 10 8 , 2.5 x 10 8 , 8.1 x 10 7 cm" 3 . We have also done 
other calculations in which we kept p ext fixed and varied T ext . These give 
very similar results. 

Fig. El shows that, right next to the star, there is a boundary layer, where 
the density rises sharply as one goes into the star and the temperature drops 
suddenly. We do not analyze this region. Once we are outside the boundary 
layer, the gas behaves very much according to the analytical solutions dis- 
cussed in Section l3~2l Starting just outside the boundary layer and extending 
over a wide range of radius, the numerical solution exhibits a self-similar be- 
havior with power-law dependences of the density, temperature and angular 
velocity. This region corresponds to the self-similar solution of MN01. There 
are, in fact, two z ones, an inner two-temperatu re zone, and an outer one- 
temperature zone (|Medvedev fc NaravanlfeoOlah [see also Eq. {ST}]. The 
most notable feature of this region is that the density, temperature and an- 
gular velocity of the numerical solutions are completely independent of the 
outer temperature and density, as predicted by the analytical solution. The 
slopes of the numerical curves also agree well with the analytical scalings. 

At a radius -R ma tch ~ 5 x 10 4 ...2 x 10 5 R g [depending on the outer pressure, 
see Eq. (|39~j) ]. solution 1 merges with solution 2 [Eq. described in the 

previous section. Here, the solution does depend on the outer boundary 
conditions, and it scales roughly according to the slopes derived analytically. 
At even larger radii R > R ext ~ 3 x 10 5 ...2 x 10 6 R g [see Eq. (gDJ)], the flow 
matches onto the ambient external medium. In this region we have solution 
3 [Eq. (jJTJ)] described in the present section. As expected, out here only 
the angular velocity and the radial velocity vary with radius. Both have the 
scalings predicted for solution 3. 
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Figure 3: Profiles of density (top left panel), temperature (top right panel, 
the electron temperature is the lower curve on the left and the proton tem- 
perature is the higher curve), angular velocity (bottom left panel), and radial 
velocity (bottom right panel), for four numerical solutions of the full height- 
integrated differential equations. The four solutions correspond to different 
values of the density of the external medium: p ext = 2.5 x 10 9 , 8.1 x 10 8 , 2.5 x 
10 s , 8.1 x 10 7 cm -3 . The first and fourth solutions are labeled 1 and 4, re- 
spectively. The temperature of the external medium and the accretion rate 
are kept fixed in all the solutions: T p ext = T e ext = 10 8 K, m = 2 x 1CT 5 . The 
analytical slopes of the three self-similar solutions are shown for comparison 
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6 Boundary layer solution 



Finally we discuss the last piece: the structure of the boundary layer. Un- 
like the previous cases, the solution to the boundary layer (BL) cannot be 
obtained in a self-similar form in terms of the radial coordinate R. The 
structure of the BL is intrinsically non-self-similar in R since all the gas pa- 
rameters (e.g., the temperature, gas density, etc.) change very dramatically 
over a relatively short radial region: R* < R < 2R* , as can be seen from 
the numerical solutions discussed in the previous sections. For instance, the 
density nearly diverges as one gets close to the star surface whereas the tem- 
perature decreases to the values well below the virial temperature. Such a 
behavior, however, suggests to look for a self-similar solution in terms of 
the distance from the stellar surface, i.e., in terms of D = R — R*. In this 
calculation we neglect the effects of radiation transfer and, especially, the 
Comptonization. They may be important in hot regions, but will unlikely 
strongly affect the flow closer to the star, where the temperature of the gas 
falls below fewxlO 9 K or so (see discussion in Section 17. 2. 3|) . 

Unlike all previous cases, here we cannot neglect the radial (infall) veloc- 
ity. For simplicity, we consider the one-temperature case. The generalization 
to the two-temperature case is straightforward: it follows from the equality of 
the Coulomb energy transfer rate and the heating/cooling rates (see Section 
13. 2j) and will be discussed in a separate paper. We again use the height- 
integrated hydrodynamic equations, but now written in the approximation 
that R = R* + D with D < R* : 



—M = 4nR*pv, 



(43) 



(44) 



4™ (pel) ^^ = J~ MQRl 



(45) 




(46) 



where 




(47) 



(48) 
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and in equation (pfHJ) we took into account that the radial velocity is negative 
(inward). 

As we mentioned above, we are looking in the solution which is self- 
similar in D, that is the temperature, density, angular and radial velocities 
are expressed as power-laws. In addition, we must satisfy the boundary 
condition at the star surface: f2 = fi*. Thus, we readily conclude that 
Q oc D° (otherwise it is either zero or diverges at D = 0, i.e., R = R*). 

Let us now consider equation (jUJ). First, we note that the rotation is 
sub-Keplerian, fll -C Vt 2 Kif so that we neglect this term. Next, we cast it into 
the form: 



M^HG^t'^ (49) 

where vg^ = V2£lk*R* is the free-fall velocity that near the stellar surface. 
We now make the following assumptions, which consistency with the obtained 
solution must be checked a posteriori: (i) the flow is always subsonic, v 2 < c^, 
and (ii) c 2 s grows with D slower than linearly (for c 2 oc D 2/3 we should have 
2(3 < 1), then the second term in the second brackets is sub-dominant and 
may be neglected as well. With these assumptions, the equation simplifies 
to 

pel = p = const., (50) 

that is, the pressure is constant throughout the boundary layer. 

Considering equation (jlSJ, we notice that since Q =const., its radial 
derivative vanishes, dVt/dD = 0, and the equation simply defines the angular 
momentum flux, J = MVt^R 2 . By the same token, the heating rate in equa- 
tion (J46|) vanishes. Together with the continuity equation (|43|) . the energy 
equation reads, 

M ( 1 ale 2 , c\ dp\ 2 

— 77; = Qs,nrP c s (51) 



AnR 2 V7 - 1 dD p dD / 

The system of equations (}4~Hj) . (joTIj) . (JSHj) . together with Q = admits 
the following self-similar solution: 

p = p Q d- 2/ \ T = T d 2/5 , v = v d 2/5 , Q = Q d°, (52) 

where we used the dimensionless distance d = D/R*. The constant factors 
are (recall that c 2 = kT/m. 



pi 
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where we denoted 

B = 107rQ ff ^(^) 5 (54) 

and p ou t is the pressure on the outside of the boundary layer, where it should 
match to the pressure on the inside of the hot brake flow. Using equations 
of section 13.21 we can calculate this pressure to be equal to 

3 , 2 , Qff,NRC 5 f RgV f m e \ 1/2 

With the solution we check that the assumptions made in order to simplify 
equation (|4T)j) are consistent: indeed v 2 /c 2 oc d 2 ^ 5 — > as d — > and 2/3 = 
2/5 < 1. There is a little subtlety with Eq. (|49p. namely the condition 
pc 2 =const. ensures that the two leading terms cancel each other exactly and 
it looks like we need to keep higher order terms, e.g., those that have been 
dropped out. Apparently, it is important to understand that the presented 
self-similar solution is approximate: the flow is pressure-dominated but with 
pc 2 ^constant. A small deviation from the exact equality is necessary to 
compensate for the remaining next order terms. 

Using our computer code, we calculated the structure of the bound- 
ary layer numerically and presented it in Figure EJ Note the remarkable 
agreement of this numerical solution with the theoretical one: Q =const., 
p «const., and p, T p , v follow the predicted scalings. 



7 Properties of the hot brake accretion flow 
7.1 Stability of the flow 

The existence of a mathematical solution to hydrodynamic equations does not 
imply that the corresponding accretion flow may realize in nature. If the flow 
is unstable to a certain type of instability, this instability may dramatically 
change the structure of the entire flow, or even prohibit it from being realized 
in nature. Here we perform the linear stability analysis of the hot brake flow. 



7.1.1 Stability to winds/outflows 

It is known that the Bernoull i parameter of the accreting gas in BH ADAFs is 

positive for a wide range ofr (|Naravan fc Yill994l l ll995atlNaravan. Kato fc Honmal 
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Figure 4: The structure of the boundary layer. Log of p, T e , T p , f2, v and 
p = pc 2 s as functions of Log(D/R*) are shown for rh — 0.01, a = 0.1, s = 0.1. 



24 



1997|), and it has been suggested that the positi ve Bernoulli parame t er may 



trigger strong winds or jets in the se systems (jNaravan fc Yilll994 Il995al 



Blandford fc Begelmanlll999L but see lAbramowicz"^t""aI"l2000[ ). flgumenshchev Sz Abramowicd 



1999, 2000) confirmed with numerical simulations that strong outflows are 
produced from BH ADAFs when a ~ 1. 

Normalizing the Bernoulli parameter, Be, by (f2^i?) 2 , and using equa- 
tions (|2 lj) . (J24|) . we find that the self-similar hot brake flow has 



Be 1 I 1 o , l n 2n2 n 2 D 2 , 7 2 

age = <U + 2 ni - + — c * 

2 7 -3 



3( 7 -1) ! 



(56) 



where 7 is the mean adiabatic index of the flow. 

The right hand side of equation (jBTIjl can be either positive or negative, 
depending on the value of 7. Hence, we find that the gas is gravitationally 
bound and unable to flow out in a wind (i.e. b < 0) if the adiabatic index 
satisfies 

7 > I (57) 

that is, the accretion flow can produce a wind and/or a collimated outflow 
only if 7 is well below that of an ideal has (7 = 5/3) and is stable to such 
outflows if 7 > 1.5. Normally, we expect 7 to be close to 5/3 for the accreting 
gas. 



7.1.2 Stability to convection 

It is well known that if the entropy increases inward in a gravitationally- 
bound non-rotating system, the gas is convectively unstable; otherwise the 
flow is stable. The specific entropy profile in the settling accretion flow 
around a NS can be readily calculated from equations (JU),© using (|2*T|) . (PH) . 
This gives 

ds_ _k 1 d_ L / jj \ _ k 2 7 - 3 1 

dR ~ m p 7 - ldR H VP 7 ™ 1 / m v 7 ~ 1 R 

We see that the entropy increases outward for 7 > 1.5 and inward for 7 < 1.5. 
Hence if 7 > 1.5 the flow is stable against convection, while if 7 < 1.5 the 
flow is convectively unstable. 
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In the presence o f rot ation, the analysis is a littl e more complicated. 
Naravan et al. ( 2000f ) and Quataert fc Gruzinov ( 2000f ) discuss the general- 
ization of the Schwarzchild criterion for accretion flows with rotation. If the 
gas motions are restricted to the equatorial plane of a height-integrated flow, 
convective stability requires the following effective frequency to be positive: 



N 2 + k 2 > 0, 



(59) 



where N is the Brunt- Vaisala frequency and k is the epicyclic frequency, 
k = f2 for n oc i?~ 3 / 2 . For a power-law flow with p oc R~ a and Q(R) = 
sttx oc -R~ 3 / 2 with s 2 = 1 — (1 + a)cg, this criterion may be written as follows 
Naravan et allEoonl more discussion^ 



sec 



N 2 



o 2 



-[( 7 +l)-a( T -l)] 



[1 + <>■)<% 

7 



+ 1 > 0. 



(60) 



Since for the self-similar settling solution a = 2, the stability criterion (jSOjl 
becomes 



A'. 



eff 



O 2 

7 



(27 - 3) > 



(61) 



which yields that the flow is convectively stable if 



(62) 



This condition is different from the stability criterion against outflows, given 
in equation (|57j). 

Fo llowing; the techniques developed bv lQuataert Sz Gruzinov (2000). iNaravan et al 
(2000) have also presented a more general analysis of convection in a self- 
similar accretion flow. This analysis, which does not restrict motions to 
lie in the equatorial plane, assumes that and c s are independent of the 
polar angle 9 (as is valid f or a marginally conv ectively stable system, cf 
Quataert fc Gruzinov 2000). iNaravanet alJ » find that the most unsta- 
ble region of the flow is near the rotation axis, 9 — 0, ir. They show that the 
marginal stability criterion for this polar fluid coincides with the condition 
for the positivity of the Bernoulli parameter. That is, a flow which is convec- 
tively stable at all 9 has a negative Bernoulli parameter, while a flow which 
is convectively unstable for at least some values of 9 has a positive Bernoulli 
parameter. (The Bernoulli parameter itself is independent of 9.) 
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We have verified this result for the solutions presented in this paper. 
Specifically, when we apply to our solutio n the more gen e ral co nvective sta- 
bility criterion given by equation (A9) of Naravan et al. ( 2000l ). we recover 
the condition (|57jl above. 



7.1.3 Thermal stability 

Not all hot accretion flows are stable. For instance, the cooling-dominated 

SLE solution has been show n to be thermally unstable ( Pir anil 19 781 : Wandel &: Liang 
19911: iNaravan fc Yilll995bh and, hence, unlikely to exist in nature. More gen- 
erally, it has been shown that any accretion flow in which heating balances 
cooling is thermally unstable if the cool ing is due to bremsstrahlung emission 
( Shakura fc Sunvaev 19761 : Piran 1978f ). The ADAF solution, on the other 
hand , is kn own to be thermally stable ( Naravan fc Yi 1995bl : iKato et"aL 



19961 1997 ). In this solution, cooling is weak (ideally zero), and so the 
thermal energy of the flow is not radiated but is advected with the gas 
(hence the name). The CDAF is also believed to be stable, since in this 
flow again the thermal energy is advected by convective eddies and is either 
carried into the black hole or is r adiated near the outer boundary of the flow 
(Ba ll. Naravan &: Quataertll200ll ). In contrast to these radiatively inefficient 
flow, the hot brake flow is cooling-dominated. Energetically, this flow is very 
similar to the SLE solution, since the heat energy produced by viscous dis- 
sipation is radiated locally via bremsstrahlung. Therefore, in analogy with 
the SLE solution, one might expect the flow to be thermally unstable. How- 
ever, this is not necessarily the case, as we show in this section (see also, 
Medvedev fc NaravanlEnOlbl ). 

The physics of the thermal instability is simple (Field 19651 ). Suppose a 
system is in thermal equilibrium, so that the rates of heating and cooling 
per unit volume are equal: Q + = Q~ . For simplicity let us take the heating 
and cooling rates to be functions of only the local temperature: Q + oc T a , 
Q~ oc T 13 (a, (3 > for concreteness) . 

Suppose, with increasing temperature, the cooling rate rises faster than 
the heating rate, i.e., (3 > a. Then a local perturbation which causes a small 
increase in the temperature will result in a net cooling of the gas: Q~ > Q + . 
This will cause the temperature to return to its equilibrium value, which 
means that the gas will be thermally stable. (It is easily seen that this 
is true also for a small decrease in the temperature.) On the other hand, 
if a. > [3, the gas is thermally unstable. For instance, if the temperature 
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decreases slightly, cooling becomes stronger than heating and the system 
deviates from its equilibrium in a run-away manner. 

When a gas is hot, one should always include the effect of thermal conduc- 
tion, — being usually a strong function of temperature, thermal conduction 
may (and will) affect the accretion flow structure and, especially, its stability. 

It is easy to see that thermal conduction will tend to reduce the thermal 
instability. An unstable thermal mode of wave-vector k consists of a growing 
temperature perturbation of wave-length 2ir/k. Thermal conduction tends to 
smooth out this temperature perturbation through heat diffusion. If the rate 
at which the temperature perturbation grows is smaller than the rate at which 
it is smoothed out by conduction, then the instability will be suppressed and 
the mode will be stable. Otherwise, the mode will continue to grow, but at 
a somewhat reduced rate. 

The rate at which fluctuations are smoothed out by conduction depends 
on the spatial scale of the perturbation. The smaller the scale (i.e. the larger 
the value of k), the faster the conduction, and the greater the stabilizing 
effect. Thus, we expect conduction to stabilize thermal modes with k greater 
than some critical k cvlt . Our task in this section is to estimate k cvlt through 
a quantitative analysis. If we find that k cri tR 3> 1, then we conclude that 
the flow is thermally unstable. On the other hand, if we find that k CT1 tR < 1, 
we may reasonably claim that the flow is thermally stable. Technically, for 
k ~ we need to carry out a global analysis rather than the local analysis 
presented in this paper, but this is beyond the scope of the present section. 

Let us write the heat flux q due to thermal conduction as 

<?cond = -kVT, (63) 
where k is the thermal conductivity c oefficien t . The rmal conductivity in a 



dense, fully ionized gas is given by the S pitzerl (J1962J) formula, 



n Sp « 1.3nk B v T X ^ 6.2 x l(T 7 T e 5/2 erg/(s K cm). (64) 

Here vt = {kBT e /m e ) l l 2 is the electron thermal speed, T e is the electron 
temperature (T e = T for a one-temperature plasma), ks is the Boltzmann 
constant, and 

A ~ 10 4 T e 2 /n cm (65) 

is the electron mean free path. Note that A is independent of the mass of the 
particle. 
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In the collisionless regime, i.e., when the mean free path of an elec- 
tron becomes comparable to or larger than the temperature gradient scale 
A > T e /|VT e |, equation (JUHj) for the heat flux is no longer valid. For 
an unmagnetized plasm a, the heat flux takes the following saturated form 
f)Cowie fc McKeelll977h . 



g sat ~ -Cpcl sgn(VT), (66) 

where C ~ 5 is a numerical constant whose exact value depends on the 
particle distribution function. This result is not relevant for our problem 
since our plasma is magnetized. 

For a collisionless magnetized plasma, thermal conduction is anisotropic. 
Electrons stream freely along the field lines, and the parallel heat flux re- 
mains the same as for the unmagnetized case described above. However, 
the transverse heat flux is greatly reduced because electrons are tied to the 
field lines on the scale of the Larmor orbit. In fact, if the field is uniform 
and homogeneous, the perpendicular thermal flux is identically equal to zero 
since electrons cannot move across the field lines. In a tangled field, however, 
electrons can jump from one field line to another and thus conduct heat per- 
pendicular to the field. Since we are dealing with a turbulent accretion flow 
with a tangled magnetic field, this is the regime of interest to us. 



The physics of this regime of conduction has been discus s ed bv|Rechester fc Rosenbluth 
( 19781 ) : IChandran fc Cowlevi (|l998h : lMedvedev fc Naravanl (|2001bh . whoide~ 



tified two important effects, which we discuss now. 

First, since particles can move freely only along field lines, the character- 
istic effective mean free path is set by the correlation scale of the magnetic 
field Ib- In a hot accretion flow this scale is not known in general. However, 
it is likely that turbulent motions in the flow occur on a scale comparable to 
the local radius R, since this is the only characteristic scale in the problem. 
Very likely, the turbulent magnetic field will also have the same scale Ib ~ R- 
We parameterize this scale as Ib = £,R- We expect £ < 1 because turbulent 
fluctuations cannot have a scale larger than the local radius of the flow. We 
assume £ ~ 0.1 throughout the paper. 

Second, the magnetic field is inhomogeneous. Therefore, only a fraction 
d < 1 of the particles will be able to pass though the magnetic mirrors that 
will be present in the field, and it is only these particles that transport energy 
beyond a distance ~ Ib- For magnetic field strength fluctuations 5B ~ (B), 
the fraction of free streaming particles is estimated to be $ ~ 0.3. 
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Typically, hot accretion flows are highly collisionless, i.e., A 3> R > Ib- 
Therefore, we can write the thermal conduction coefficient as 



k b ^ nkEVrls^ — 10 nk B VTR£,-i$-i, (67) 

where £_i = ^/10 _1 and = t9/10 _1 . Let us write the conductive heat 
flux in a form similar to that used for the viscous stress, namely 

c 2 dc 2 

<?cond = -^c-^-p-T-, (68) 

where the dimensionless coefficient a c is analogous to the Shakura-Sunyaev 
viscosity parameter a, and is given by 

a c ~ ^^-10~ 2 ^-i. (69) 

Here we have used the fact that vt — c se and H/R ~ c s /vq ~ c s /QkR, 
where H is the accretion disk scale height (in hot flows, H ~ R) and Vg is 
the free-fall speed. 

To study the thermal stability of an accretion flow, we need to include ad- 
ditional physics, nam ely the effects of shear and rotation. We use the shearing 
sheet approximation dGoldreich fc Lvnden- Bellll96,4ljulian fc ToomrJl96f 



Goldreich fc Tremainelll978l ). which is a convenient way of introducing the 



relevant physics without unnecessary technical complications. Furthermore, 
we carry out a local WKB analysis under the assumption that the wavelength 
of the perturbation is much smaller than the radius. 

The shearing sheet model approximates the flow as locally flat, neglecting 
the effects of the flow curvature. Conventionally, the shearing sheet coordi- 
nates are Cartesian with x, y, z corresponding to the radial, azimuthal, and 
vertical directions, respectively, centered on some point in the flow at a ra- 
dius R. These coordinates are appropriate for describing the motion of a 
parcel of gas whose geometrical size is small compared to the local radius, R, 
of the flow (i.e., x, y, z R), so that the effects of geometry and curvature 
are insignificant. We neglect viscosity in the azimuthal momentum equation; 
of course, we do include viscous dissipation in the energy equation, where it 
plays an important role. 

It is convenient to compare the wave-vector k of a perturbation with 
1/R and the frequency of a mode with the local Keplerian frequency Qk — 
^/GM*/R 3 , where M is the mass of the central object. The shearing sheet 
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approximation is accurate for "local" small-scale perturbations with kR ^> 1. 
Perturbations with kR ~ 1 are global; their properties may be understood 
only through a global stability analysis, usually numerical, which we do not 
attempt here. 

We consider a shearing gas flow with unperturbed velocity given by 

V (x) = 2Axy, (70) 

where 2 A = dVo/dx is the shear frequency and "hat" denotes a unit vector. 
Note that we have neglected the radial velocity in the equilibrium flow since 
this component of the velocity is significantly smaller than the azimuthal 
velocity. To include the effect of rotation we assume that there is a Coriolis 
acceleration, described by an angular rotation frequency fl = Q z. The 
vorticity and epicyclic frequency are then given by 

2B = 2A + 2tt, K epi = 2(OB) 1/2 . (71) 

The hot brake flow satisfies the Keplerian scaling, f2 oc i?~ 3 / 2 . Therefore, we 
have 2A = -(3/2)Q, 2B = Vt/2 and K cpi = Q. 

We assume that perturbations in the flow have structure only in the x di- 
rection, and we ignore motions in the z direction. We write the perturbations 
(represented by primes) in the velocity, density and sound speed as 

V'(x,t) = u(x,t)x + v(x,t)y, (72) 
p'(x,t) = p a(x,t), (73) 
cffof) = a 2 (x,t), (74) 

where po an d c 2 are the equilibrium values of the density and the square of 
the sound speed. Note that we define c s to be the isothermal sound speed, 
so that the pressure is written as p = pcj.. By considering perturbations 
of the basic hydrodynamic equations, namely the continuity, radial momen- 
tum, azimuthal momentum and entropy equations, we obtain the following 
four linearized equations (note that in the presence of conduction the energy 
equation has an additional contribution from the divergence of g CO nd, and in 
equilibrium, the total heating is equal to the total cooling), 

a + &=°- (75) 

du ^ 9 da da 2 , 

w- 2f! " + % + &=°- (76) 
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dv 
dt 



+ 2Bu = 0, 



Po da 2 



>do 



7 -i at ~ poCs m 



(77) 
(78) 



fix dx 2 
d/dt since the inflow velocity 



where we have used d/dt = d/dt + V 0x d/dx 
Vox is set to zero in our approximation. 

For the heating and cooling rates, we make use of "realistic" expressions 
that represent the physics of viscous accretion flows. Thus we write 



or 



a 



P^_ 



dVr 



Oy 



dx 



4aA 



2P C s 



n 



Q- 



2\ n 



K 



-CP 2 (cl) 



(79) 



where a ~ 0.1 is the standard Shakura-Sunyaev viscosity parameter, Vq v is 
the ^/-component of the unperturbed velocity, and C is a constant. We leave 
the index n in the cooling function unspecified for now, but we note that 
n = 1/2 corresponds to non-relativistic free-free (bremsstrahlung) cooling. 

We assume that the perturbations in equations (J75|) are of the form 
exp(— iut + ikx). Substituting in the above equations and solvi ng;, we obtain 
the following dispersion relation (jMedvedev fc NaravarJl2l)01bh : 



CO 



7 



+ 



i\n 



1) ik 2 K 
— + - 
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Tcool 



0J + 



''"cond 

i(2B/A 



Tcool 



2„2 



k 2 c 



where 



^cool 



n 



K 



AA 2 a 



9as 2 n « 



is the cooling (heating) time of the gas and 

Q K R 2 /a c c 2 



Tcond 



(80) 



(81) 



(82) 



is the conductive time scale. 

The dispersion relation ()80|) corresponds to purely radial perturbations. 
The same relation can be used also for perturbations in the vertical direc- 
tion, except that we must set k 



0. Perturbations in the azimuthal 
direction are more complicated. Because of the shear, a non-axisymmetric 



cpi 
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wave packet is distorted as a function of time, and must be analyzed by spe- 
cial techniques which are beyond the scope of this study (see, e.g., Toomrel 
19771: lOoldreich fc Tremainelll978h . 



Equation ()80|) is a fourth-order polynomial and has four roots correspond- 
ing to four modes. A flow is unstable if any of the four modes grows with 
time, i.e. if the corresponding root has Im uj > 0. One of the roots of the 
dispersion relation is always uj = 0. This root corresponds to the viscous 
mode, which in the present case is particularly simple because we neglected 
viscosity in the momentum equation. It is easy to show that if we introduce 
viscosity into the momentum equation the viscous mode would become sta- 
ble, i.e., we will obtain Im uj < 0. We do not consider the viscous mode 
further. 

Let us now neglect thermal conduction for a moment, r con d — > oo. Then 
the physics of the remaining three modes may be understood by considering 
equation (JHUj) in various limits. Consider first the limit k — > 0. In this limit, 
two of the roots are given by uj = ±K cpi , corresponding to simple epicyclic 
oscillations. In the opposite limit k — > oo, the same roots are given by 
uj = ±7 1//2 c s fc, which shows that they correspond to sound waves. In the 
absence of heating and cooling (i.e. r coo i — > oo), we can obtain an exact 
solution for these roots which is valid for all k: 

uj 2 = K 2 cpi + 1C 2 s k 2 . (83) 

This is the standard dispersion relation for sound waves in a differentially 
rotating flow. The presence of 7 is because the relevant sound speed is the 
adiabatic sound speed, 7 1//2 c s (recall that c s is defined to be the isothermal 
sound speed). 

The final root of the dispersion relation (jHUj) corresponds to the thermal 
mode. In the limit k — > 0, we obtain 

UJ = l { ri -l)^—A, ( 8 4) 

7"cool 

We see that the mode is stable (for 7 > 1) if n > 1 and unstable if n < 1. In 
the opposite limit k — > 00, we find 

■ (7-1) (2-n-2B/A) 
uj = 1 . (85) 

7 ''"cool 

Therefore, the mode is stable if n > 2(1 — B/A), i.e. n > 8/3 for our problem, 
and unstable otherwise. Note that an accretion flow that is cooled by free-free 
emission (n = 1/2) is unstable in both limits. 
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Let us now turn the conduction back on. Figure shows that the imag- 
inary parts of all three modes decrease rapidly with increasing k. For the 
particular parameters we have selected, the growth rate of the unstable ther- 
mal mode (curve 3) goes to zero at k CTit R ~ 1.5, and the mode is stable for 
all k ^> k C Y\t' 

We may also analyze equation (J8(J|) analytically. In the large-/c limit, the 
root corresponding to the thermal mode is equal to 1 

7 V 7~cool "^cond / 

Clearly, for large k, conduction stabilizes the thermal mode, for the reasons 
explained at the beginning of this section. Using the above relation, we can 
estimate the critical k cvlt above which all k are stable: 

5 2 7"cond I r\ c, B \ 13 T conc j 



Tcooi V AJ b r coo i 

where we have substituted n — 1/2 (free-free cooling) and B/A = —1/3 
(Keplerian scaling). 

We should comment here that in the theory described above, the conduc- 
tivity Kb is a quantity averaged over many field correlation lengths. There- 
fore, the results of the stability analysis are valid only for perturbations 
on scales much larger than Ib- (This is somewhat inconsistent since we 
have assumed that Ib = £,R ~ 0.1-R.) For small-scale perturbations with 
k ^> (£i?) _1 , the local magnetic field is nearly homogeneous. Therefore, 
thermal conductivity is anisotropic; its perpe ndicular component is of orde r 
of Kb ~ n/c B t> r / s ^ (for more discussion, see, Medvedev fc Naravan 2001bj ). 



while the conductivity along the field lines is much larger: 
9H _ CpcUgn(VT) Cpcl 5 



>1, 



q± -a c {pc 2 JVt K ){ dc 2 Jdx) a c pc 3 s (H/R) a c 

1 In deriving this equation from (|80J) we have used the fact that the acoustic time-scale 
is, in general, shorter than the time-scale of the thermal mode, i.e., w <^ kc s , and we 
have neglected K ep i as before. In this case we can neglect ui 2 in the second brackets, 
so that equation (|86|l readily follows. It may seem that this procedure fails when the 
acoustic and thermal time-scales are comparable. This may happen when kR ~ 1 (for 
such perturbations, the sound crossing time is of order the dynamical time) and when 
t coo i and T con( j are also comparable to the dynamical time, which is ~ ^j^- 1 - Nevertheless, 
even in this case, equation (|86|) works fairly well near the stability threshold. Indeed, at 
the threshold itself, Im uj — 0. Since further the thermal mode frequency has no real part, 
we have u> ~ near the threshold; so we may safely neglect uj 2 compared to k 2 c 2 s in the 
second brackets of equation Ij86(l . 
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Figure 5: Real part (left panel) and imaginary part (right panel) of the 
frequencies of three modes. The curves labeled 1 and 2 refer to the two 
acoustic modes, and the curves labeled 3 refer to the thermal mode. The 
dashed curves correspond to the dispersion relation (JRHj) . which does not 
include thermal conduction r cond — > oo. The following parameter values 
were used: n epi /Q K = s = \/ol), t coo \ :S = lOttj}, c s = Vt K R, 7 = 5/3. 
The mode frequencies are normalized by the Keplerian frequency. The solid 
curves correspond to the dispersion relation (fKUJ) . which includes thermal 
conduction. Here, r con d = Tcooi^, and the other parameters are the same as 
before. 
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as follows from equations (|68|) and (|66|) . Thermal instability along the mag- 
netic field lines is then suppressed more strongly than in the analysis pre- 
sented above. 



Thermal conduction in the hot brake flow is enormous (IMedvedev fc Naravan 



2001bh due to the very high temperature and large mean free path of parti- 



cles. The estimated value of thermal conduction is 

a c ~ few x l(T 2 f_itf_i (88) 

Is this level of thermal conduction enough to stabilize the thermal mode? 

The relevant stability criterion is given in equation (J87j) . However, before 
we apply this criterion, we need to allow for the fact that thermal conduc- 
tion in the flow is so strong that it modifies even the equilibrium structure 
of the flow. In particular, th e cooling time (JHTj) is modified and becomes 
()Medvedev fc NaravanlboOlbh 



^^"'4"'' (89) 

where we have assumed that a c ^> as 2 , which is reasonable for typical pa- 
rameters, e.g. a ~ 0.1, s ~ 0.1, a c ~ 0.1. From equation (J%2~)) and using 
the self-similar solution from Section 13.21 we obtain the thermal conductive 
time 2 

T omd = — nj c 1 , (90) 

OL c . 



Substituting r coo i and r cond into the stability criterion (JHTj) . we then find 



9(9as 2 + 2a c 



1/2 , 13 




- ~ 1.2, (91) 



that is, thermal modes with kR > 1 are stable. This suggests that the hot 
brake flow with thermal conduction is stable to the thermal instability. 

Whether the mode kR = 1 itself is stable or not cannot be reliably de- 
termined from our local analysis. A global stability analysis is necessary to 
properly account for the effects of geometry and curvature, but this is beyond 
the scope of the present analysis. 



2 Alternatively, we recall that in the hot settling flow H/R ~ 1 and H/R ~ c s /vg 
c s / (V2QkR)- Then r con( j readily follows from equation l|82() . 
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7.2 Observable properties 

7.2.1 Spin-Up/Spin-Down of the Neutron Star 

The rate of spin-up of the accreting NS is given by 

4 = J - MM(R*)Rl ^ -43s 3 a 2 M Edd n K (R*)Ri (92) 
at 

where I* is the moment of inertia of the NS. We have made use of the fact 
that | J| \MVL(R*)R 2 \ for the self-similar solution, and used equation (fTUj) 
for J. The negative sign in the final expression implies that the accretion 
flow spins down the star. The above equation is for an unmagnetized NS. If 
the NS has a magnetosphere, the inner edge of the accretion flow is at the 
magnetospheric radius, R m . In this case, let us define s by = sflxiRm) = 
sf2x(-R*)(-Rm/-R*)~ 3//2 - Substituting this in equation (J9~2~j) with = constant 
and integrating, we obtain 

j / p \ -3/2 



y/T+t/r' 86s 2 a*M Edd R*\R 

where so = s(t = 0). The same result is valid for an unmagnetized NS by set- 
ting R m = R*. The quantity r is the characteristic spin-down time of the NS. 
For a spherical NS of constant density, 7* = 2M*R 2 /'o = (0.8 x 10 33 g)mR 2 . 
Substituting this expression, we obtain the spin-down rate P*/P* = r _1 with 

(P \ -3/2 / p \ -3/2 

^r) s = 2xlO^ M (f 1 yr. (94) 

Note the remarkable fact that the spin-down time scale is independent of the 
mass of the NS, and the mass accretion rate! For the magnetic case, the rate 
depends on the radius ratio R m /R*. 

It is customary to express the spin-down rate as P*/P 2 . Writing 



2vr (R 



3/2 

p * = s n K (R«) lx J (95) 

and Qk(R*) — 10 4 m7i rad/s, where R* = 3R 9 and m 14 = M*/(1.4M Sun ), 
we obtain 

p 

~ 2.4 x l(r 10 m->V s~ 2 = 2.7 x Kr 12 mr> 2 3 s 3 5 s -2 , (96) 
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where s . 5 = s/0.5. This spin-down rate is in good a greement with observa- 



tional data on the spin-down of X-ray pulsars for which I Yi. Wheeler fc Vishniac 



( 19971) invoked ADAFs: 4U 1626-67 has P/P 2 « 8 x 10" 13 s~ 2 and P = 7.7 s; 
OAO 1657-415 has P/P 2 « 2 x 10~ 12 s~ 2 and P = 38 s, and GX 1+4 has 
P/P 2 « 3.7 x 10~ 12 s -2 and P = 122 s. Since the spin-down rate is quite 
sensitive to a and s, the observed data in individual systems can be fitted 
by small adjustment of these parameters. 

7.2.2 Luminosity 

In computing the luminosity of the accretion flow, we must allow for the 
energy release in both the boundary layer and the self-similar settling zone. 
We calculate their luminosities separately. 

Radiation from the self -similar hot brake flow may be calculated following 
the methods described by Popham fc Naravan ( 1995^ 1 for a thin disk. This 



method assumes that the luminosity at a given radius is determined by the 
local viscous energy production. This is a legitimate approximation for the 
settling flow in which q~ = q + . Keeping only the dominant terms, we find 

r ^s = ^— (l~js) + M / — , (97) 



Rin JP in P 

where R in = R* + A BL is the inner radius of the self-similar zone and A BL <^ 
P* is the thickness of the boundary layer. Here the first term inside the 
parentheses represents the luminosity associated with the potential energy of 
the infalling gas, L pot , the second term — js is the luminosity associated with 
the rotational energy extracted from the star, L TOt (note, j < in the self- 
similar solution), and the final integral is the "enthalpy correction", L ent h. 
Using the analytical solution (J2l"j) - (f24j) and assuming p out = for simplicity, 
we obtain 

Try 

L P ot = M Edd c 2 — , (98) 

2 

L TOt = 43M Edd c 2 ^ 4 , (99) 
L enth = -M Edd c 2 — . (100) 
Note that the leading terms in L pot and L enth cancel each other exactly. The 
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luminosity of the self-similar flow is thus 

L ss ~ 6.2 x 10 34 mr3 1 m_ 2 sg -1 + 8.9 x 10 33 mr^a 2 01 s 4 01 ergs" 1 , (101) 

where rh_ 2 = m/0.01, So.i = s/0.1, r 3 = r*/3,and we have assumed s<l. 
Note that luminosity not associated with dissipation of rotational energy, 
represented by the first term in equation (jlUlj) . is much less than the com- 
monly assumed ~ GM*M / R*. This is because the negative enthalpy term 
has a large magnitude, as a result of the fact that the settling flow is akin to 
a pressure supported, quasi-stationary atmosphere. 

The second term in equation (jlUlj) is the luminosity of the settling zone. 
Since the self-similar solution for this zone is independent of m, the luminosity 
too shows no m dependence. Indeed, the luminosity remains finite even as 
rh —>■ 0. How is this possible, and where does the energy come from? The 
answer is that the luminosity of the settling zone is supplied by the central 
star. As the star spins down, it does work on the accretion flow and the 
energy released comes out as bremsstrahlung radiation. 

The boundary layer luminosity requires a different method of calculation 
since viscous energy production is negligible in this zone: Q ^constant, and 
so q + oc (dQ/dR) 2 ~ 0. As the accreting gas cools in the boundary layer, 
starting from a nearly virial temperature ~ 10 12 K on the outside down to 
the NS temperature ~ 10 7 K near the surface, the thermal energy in the 
gas is emitted as radiation. To estimate the luminosity, we use the energy 
balance equation, which is the sum of equations (jlj),©: 

pv dc 2 2 dp 7 dc 2 dP 
~ q ~ = —dR - VC °dR = — pV dR ~ V dR- (1 ° 2) 

We can neglect the dp/dR term because the pressure p is essentially con- 
stant in the boundary layer. To obtain the luminosity we integrate over the 
boundary layer 

Lbl = [ q~^R 2 dR = - f -^AtxR 2 P v < ^- dR = — - — M (Ac 2 ) . (103) 
J J J - 1 dR 7-I 

Since c 2 starts from nearly virial value and reaches close to zero, Ac 2 s ~ 
GMJR*. More precisely, Ac 2 , = c 2 A(6 p + 6 e ) ~ c 2 9 p0 r~ 1 . Therefore, the 
boundary layer luminosity is 

Lbl = -^-M EAd c 2 — w 1.7 x 10 36 mm^\ (104) 
7-1 or* 
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where we have assumed 7 = 5/3. The total luminosity of the system is 
L = L S s + L BL . Note that L BL has the same scaling as the first term in L ss 
[Eq. ([101)1 ] and dominates the latter unless s — ► 1. 

7.2.3 Effect of Comptonization 

Using the self-similar solution ([21 |) .([24 p . we may readily estimate the electron 
scattering optical depth and the ^-parameter. 3 The optical depth is 

T es ^ pK es R ^ lO^sV" 1 ~ a .l4 1 7 " 1 ' ( 105 ) 

where K es = OTjvn v is the electron scattering opacity for ionized hydrogen. 
Since r > 3, we see that r es < 1/3 for reasonable parameters and the radiation 
is optically thin to electron scattering. The ^/-parameter is 

y = l69 2 e T es ~ 2 x 10 6 as 2 r~ 2 ~ 2 x 10 3 a .i4i r ~ 2 - ( 106 ) 
The radius at which y ~ 1 is 

r c ~ 45^^0.1- (107) 

Above this radius the inverse Compton scattering is small and the self-similar 
solution is valid. For r < r c , however, Comptonization is important and 
the electron temperature profile will be modified from the self-similar form. 
Since the electron-proton collisions are relatively weak (the plasma is two- 
temperature), other quantities, e.g., the density, proton temperature, etc., are 
unaffected. Comptonization is unimportant for low-viscosity flows, a < 0.01 
around slowly rotating NSs, s < 0.01, because then r c < r*. 

7.2.4 Spectrum 

We now estimate the spectrum of radiation emitted from the hot brake ac- 
cretion flow. Let us neglect inverse Compton scattering for the moment. The 
relativistic bremsstrahlung emissivity is approximated as 

t v oc p 2 exp — {hu/kT e ) erg cm -3 s -1 Hz -1 . (108) 

3 Here we just estimate where the effect of Compt onization becomes significant. Fo r 
better analytical approximati ons see, for instance, iDermer. Liang fc Canfieldl l|l99l|) : 
iTitarchuk &: Lvubarskiil (^)95). (In the latter paper, the expression for y is not given, 
but it can be inferred using equation [24]: y = r es [(a + 3 W(1 + 6)+ AdH a 9 2 \.) Comp- 
tonization of free- free radiation has also been considered bv lTitarchukl l|l989l) . 
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Therefore the luminosity per unit frequency is 

POO POO 

L u (x p 2 e- hu/kn 2nR 2 dR(x r 3 e~ vt dt oc u 2 T(-2, u/u m ), (109) 

where T(a,z) = f°° £ 0_1 e~*cf£ is the incomplete gamma-function and v m = 
kT e (R*)/h is the maximum frequency. Above u m the spectrum falls exponen- 
tially and below u m it is nearly flat. We may, thus, replace the exponential 
in the integral with a square function which is equal to unity for v < v m and 
for v > v m . With this approximation 

L sw^y (11 o) 



1 V m \ V, 



where L$s — J L u dv is the total luminosity of the self-similar flow, rep- 
resented by equation (jl01|) . The break frequency, v mi is roughly given by 
hv m ~ 2.7 MeV for a typical electron temperature T e max ~ 10 10 ' 5 °K [cf., 
equation (J2U)]- At a typical x-ray energy, hv ~ 3 keV, the observed luminos- 
ity per decade is 

vL v ~ 1.7 X 10 31 m<4 if ergs" 1 , (111) 

i.e., vL v ~ 1.5 x 10 32 for a 300 Hz neutron star (s .i ~ 1.6). The spectrum 
is very hard with a photon index of order unity. Therefore, the luminosity 
per decade is much greater at higher photon energies and may be as high as 
~ few x 10 34 - 10 35 erg/s at hv ~ MeV. 

As shown in the previous section, Comptonization becomes important 
below the radius r c . At r c , y « 1 and the electron temperature is 

T e (r c ) ~ 2.7 MeV/v^ ~ 400a^ /4 %i /2 keV. (112) 

For r < r c , the electron temperature will be determined self-consistently by 
Compton cooling rather than by bremsstrahlung emission. Computing the 
spectrum from this region is beyond the scope of the paper. We also do not 
attempt to calculate the spectrum of the radiation from the boundary layer. 



8 Conclusions 

In this paper we present the analytical self-similar solution describing this 
hot brake flow. The hot brake flow exists at accretion rates as low as a few 



41 



percent of Eddington or smaller. We studied its properties and showed that it 
is stable with respect to winds/outflows and convection, as well as thermally 
stable (when the effects of thermal conduction are self-consistently included). 
The flow is subsonic everywhere. The flow is cooling-dominated, it is powered 
by the rotational energy of the central accretor which is braked by viscous 
torques. A very interesting property of the flow is that, except for the inflow 
velocity, all gas properties, such as density, temperature, angular velocity, lu- 
minosity, and angular momentum flux, are independent of the mass accretion 
rate (the flow properties do depend on the NS spin). This property implies 
that the density, temperature and angular velocity of the gas at any radius 
are completely independent of the outer boundary conditions. Therefore, 
such a flow cannot be matched to a general external medium. Hence, there 
is a "transition region", represented by another self-similar solution, which 
has an extra degree of freedom which permits the hot accretion flow to match 
general outer boundary conditions. Matching the hot flow on the inside to 
general boundary conditions at the star surface occurs via another transition 
region, known as the "boundary layer" . A self-similar solution representing 
the boundary layer has also been derived. 

Although we presented the hot break solution in the context of a neu- 
tron star, a similar accretion flow will certainly occur around a rapidly 
spinning white dwar f , and even around a spinning black hole, provided the 
Blandford fc Znaiek (11977 ) mechanism is efficient enough to power the MHD 



hot b rake flow (|Medvedev fc MenoiJl2002l : IPerna. et al.ll2003[ iMedvedev fc Murray 
20021) . 
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